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Abstract 

Background: Attitudes and policy towards smoking changed over the past years in many countries including the 
Netherlands. Generally, this led to a decrease in smoking prevalence. As demonstrated in twin and family studies, 
individual differences in smoking behavior are partly influenced by genetic factors. We explore whether the current 
change in environmental conditions has influenced the genetic architecture of smoking. This would constitute 
evidence for Gene x Environment (GxE) interaction. 

Methods: Data on smoking were available from 2 cohorts of young adult twins (18-25 year) registered with the 
Netherlands Twin Register. The first cohort completed a survey in 1993-1995 (n = 2669) and the second in 2009- 
2010 (n = 2339). Prevalence and genetic architecture of smoking were compared across cohorts using structural 
equation models in MX. 

Results: Smoking prevalence decreased from 40-51% to 22-23% between 1993-1995 and 2009-2010. Genetic 
analyses, making use of the different genetic resemblance in monozygotic and dizygotic twins, showed that the 
heritability was the same in both cohorts. 

Conclusions: The change in policy and smoking attitudes that led to a decrease in prevalence of smoking did not 
change the heritability of smoking and thus no evidence was found for GxE interaction. 
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Background 

Smoking ranks high among public health problems in 
the world. In the past years worldwide media campaigns 
discourage smoking behavior and emphasize the nega- 
tive health consequences of smoking. In a large study 
among high school students in the United States, the 
prevalence for ever smoked cigarettes did not change 
from 1991 to 1999, but declined between 1999 and 2009 
[1]. A study in Switzerland showed a decrease from 
1992 to 2007 in the prevalence of ever smokers in men 
but not in women [2]. The same phenomenon was 
observed in Finnish men: the prevalence of smoking 
decreased from 1972 to 2007. Among women, smoking 
increased throughout the years until 2002 but remained 
stable between 2002 and 2007 [3]. 
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In the Netherlands smoking policies changed in the 
past years. Cigarette packs contain warnings about 
health consequences since 2002 and advertisement for 
smoking is forbidden since 2002/2003 (before 2003 it 
was already forbidden for television and radio but since 
then also in newspapers, magazines et cetera). Smoking 
is forbidden in public transportation and at work since 
2004 and it is not allowed to smoke in the catering 
industry since 2008. In addition, media campaigns dis- 
courage adolescents to take up smoking and the price of 
cigarettes increased greatly. The prevalence of current 
smokers gradually declined from the nineties till now 
http:/ /www.stivoro.nl. 

Understanding the trends in the prevalence of cigar- 
ette smoking among adolescents and young adults 
enables policy makers to target prevention resources 
more effectively. In this paper, we compare the preva- 
lence of smoking in young adult (18-25 year old) twins 
who participated in the longitudinal survey study of the 
Netherlands Twin Register (NTR) in 1993-1995 with the 
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prevalence of smoking in young adult twins who partici- 
pated in 2009-2010. 

Twin data offer the unique opportunity to disentangle 
genetic and environmental influences. Previous studies 
based on data from the NTR have shown that smoking 
in is influenced both by shared environmental (51-56%) 
and by genetic factors (36-44%) [4-6]. The estimates for 
the importance of those factors are comparable with 
other twin studies worldwide [7-10]. A meta-analysis 
estimated that, for men and women, 37% and 55% of 
the variance in smoking initiation was due to genetic 
factors, and 49% and 24% due to shared environmental 
factors [11]. 

In this paper we look at differences in heritability esti- 
mates between two birth cohorts. Such differences could 
point to Gene x Environment (GxE) interaction, which 
is defined as a difference in the influence of genetic fac- 
tors conditional on environmental exposure [12]. The 
importance of the social context may influence heritabil- 
ity in several ways. According to previous studies, three 
complementary mechanisms can help explain GxE 
interaction mechanisms [13,14]. Under a social control 
mechanism the social forces wash out the effects of 
genetic factors in tobacco use. Norms and other social 
forces restrict variability in phenotype of genetically 
diverse people. For example, genetic influences on sub- 
stances use (smoking & alcohol) were low or nonexis- 
tent among those who were raised with a strong 
religious upbringing [15,16]. According to the social 
trigger model, genetic factors differentiate between indi- 
viduals only in the presence of social pressures. For 
example, the genetic influences on daily smoking were 
significantly higher in high school students who 
attended schools in which the most popular students 
smoked the most. In this case, the pro-smoking norms 
serve as trigger for the genetic influence. In the first two 
models, shared behavioral expectations and correspond- 
ing sanctions cause differences in heritability by block- 
ing or enabling genes to be expressed. In contrast, the 
social push model suggests that changes in social norms 
can affect the relevance of genetic influences by mini- 
mizing or maximizing "noise" that has the potential to 
overwhelm or hide the influences (non-causal). As large 
numbers of people (regardless of genetic makeup) begin 
smoking, there will be a tipping point in the distribution 
of smoking environments where entree into smoking 
becomes a primarily social phenomenon; genetic vulner- 
able persons are no more likely to begin smoking than 
genetically invulnerable persons simply because of the 
predominant social popularity of smoking (overwhelm- 
ing influence of social factors). This was observed in 
Swedish women born between 1910 and 1924: the pre- 
valence of smoking was very low and heritability was 0. 
Both prevalence and heritability increased in later birth 



cohorts when smoking among women became more 
accepted [17]. 

In this paper, we explored whether a change in envir- 
onmental conditions led to a change in the relative con- 
tribution of genetic factors to smoking. This would 
constitute evidence for GxE interaction. Data on smok- 
ing were available from two cohorts of young adult 
twins (18-25 year) registered with the Netherlands Twin 
Register (NTR). The first cohort completed a question- 
naire in 1993-1995 (n = 2669) and the second cohort in 
2009-2010 (n = 2339). Genetic analyses, making use of 
the different genetic resemblance in Monozygotic (MZ) 
and Dizygotic (DZ) twins, will be used to explore herit- 
ability as a function of environmental exposure (1993-95 
versus 2009-10). 

Methods 

Sample and measures 

Smoking data were collected in a longitudinal survey 
study of the Netherlands Twin Register (NTR). The 
NTR was established in 1987 and registers twins and 
their family members who volunteer to taking part in 
research studies [18]. Since 1991, every two to three 
years surveys are mailed to adolescent and adult twins 
and their family members. Surveys focus on health, life- 
style, demographics, psychopathology and personality. 
Previous studies have shown that the non-response bias 
in the NTR sample is rather small [19,20]. All the NTR- 
work with human subjects is reviewed by the Central 
Ethics Committee on Research Involving Human Sub- 
jects of the VU University Medical Center, Amsterdam, 
an Institutional Review Board (IRB) certified by the US 
Office of Human Research Protections (IRB number 
IRB-2991 under Federal wide Assurance-3703; IRB/insti- 
tute code NTR 03-180). 

To compare cohort effects over a 15-year time inter- 
val, data collected in surveys in 1993-1995 were com- 
pared to data which were collected in 2009-2010. Data 
on smoking were available for 2305 18-25 year old twins 
who participated in 1993. Additional data on smoking 
initiation were added from 364 18-25 year old twins 
who participated in 1995. This combined dataset con- 
sisted of data on 2669 twins: 415 MZ male (206 pairs), 
363 DZ male (178 pairs), 658 MZ female (326 pairs), 
462 DZ female (227 pairs) and 769 DZ opposite sex 
twins (377 pairs). 

The survey of 2009-2010 is ongoing. The available 
sample 18-25 year old twins consisted of 296 MZ male 
(105 pairs), 238 DZ male (74 pairs), 703 MZ female 
(269 pairs), 444 DZ female (149 pairs) and 649 dizygotic 
opposite sex twins (166 pairs). 

The same question on smoking initiation was asked in 
both surveys: "Did you ever smoke?" with answer cate- 
gories no, a few times to try and yes. Answers were 
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recoded into a dichotomous variable: yes versus no or a 
few times to try. 

Analyses 

Testing and genetic analyses were carried out using 
structural equation modeling in the software package 
MX [21]. A threshold model was used [22], which 
assumes an underlying (latent) liability to a categorical 
variable such as smoking. This liability is the sum of the 
effects of many genetic and environmental factors. It 
has a normal distribution with a mean of zero and unit 
variance. A threshold, based on the prevalence of smok- 
ing in the population divides the distribution into 
"affected" and "unaffected" individuals. A test for equal- 
ity of thresholds between groups amounts to testing for 
equality of prevalence between e.g. cohorts or men and 
women. The resemblances between twins for the liability 
to smoking were summarized with tetrachoric twin cor- 
relations. Based on the correlations a genetic model can 
be fitted to the data to estimate the contribution of 
additive genetic factors (A), common environment (C) 
and unique environmental (E). The total variance of the 
liability was constrained to be 1 for men and women. 
Additive genetic effects are the same for MZ twins 
because they are genetically (nearly) identical. DZ twins 
however, share only 50% of their segregating genes. 
Common environmental effects are thought to be the 
same for MZ twins and DZ twins, based on the assump- 
tion that monozygotic twins do not share more environ- 
ments with each other than DZ twins [23]. 

To investigate GxE, we tested whether the relative 
contributions of A, C and/or E were significantly dif- 
ferent in the two cohorts. In addition, quantitative and 
qualitative sex differences were assessed. To test for 
quantitative sex differences we tested if the contribu- 
tions of A, C and/or E differed between men and 
women. Finally, we assessed qualitative sex differences 
in the contribution of common environmental factors 
by specifying a free parameter in DOS twins that 
represents the correlation between the common envir- 
onment of men and women (Rcdos). If Rcdos equals 1, 
this indicates that qualitative sex differences are not 
detected. 

In the full ACE model (model 1) the influence of A, C 
and E factors was estimated separately for men and 
women and for both cohorts. The shared environmental 
correlation (Rcdos) in DZ twins of opposite sex (DOS) 
was estimated as a free parameter (in both cohorts). In 
model 2 the qualitative sex differences are tested (Rcdos 
= 1) while in model 3 to 5 the cohort differences for A, 
C and E were tested within sex. In model 6, quantitative 
sex differences were tested by equalizing the estimates 
for A, C and E for men and women. In models 7 and 8, 
A and C are dropped respectively. 



Modeling was done on the raw data, using full-infor- 
mation maximum likelihood. This gives a likelihood 
(LL) for each model that is fitted to the data. Twice the 
difference between two nested models has a chi-square 
distribution with the degrees of freedom (df) equal to 
the difference in df between the two models. A p-value 
of 0.01 was used to test if a less complex sub model 
explained the data significantly worse than the more 
complex model. 

Results 

A noteworthy decrease in smoking prevalence was 
observed from 1993-1995 to 2009-2010, in both sexes 
and in all zygosity groups (P < .001), as may be seen in 
Table 1. In the 1993-1995 cohort around 50% of the 
man and 40% of the women had initiated smoking while 
in the 2009-20 10-cohort those percentages decreased to 
around 20% in both men and women. The prevalence of 
smoking was lower in MZ twins compared to DZ twins 
in both sexes and both cohorts (p < .001). 

Tetrachoric twin correlations are shown in Table 2. 
Based on the differences in prevalence, different thresh- 
olds were retained for men and women and for mono- 
zygotic and dizygotic twins (4 different thresholds). In 
both cohorts, MZ correlations were higher than DZ cor- 
relations suggesting that genetic factors play a role in 
smoking initiation. The DZ correlations were higher 
than half the MZ correlations which suggest that shared 
environmental factors are also involved. 

Heritability estimates are shown in Table 3. In both 
cohorts estimates were very similar. 

The correlation between the common environment of 
men and women (Rcdos) could be constrained at 1 (Table 
4, model 2) indication no qualitative sex differences on 
smoking. Constraining the estimates for the contributions 
of genetic, common and unique environmental factors 



Table 1 Prevalence of smoking initiation in 18-25 year 
olds in 1993-1995 and 2009-2010 





18-25 year old twins in 
1993-1995 


18-25 year old twins in 
2009-2010 




N total 


N ever 
smoked 


% ever 
smoked 


N total 


N ever 
smoked 


% ever 
smoked 


Men: 


MZ 


415 


190 


45.8% 


301 


44 


1 4.6% 


DZ ss 


363 


184 


507% 


239 


69 


28.9% 


DZ os 


384 


218 


56.8% 


228 


62 


27.2% 


Total 


1162 


592 


50.9% 


768 


175 


22.8% 


Women: 


MZ 


658 


224 


34.0% 


728 


128 


1 7.6% 


DZ ss 


462 


195 


42.2% 


445 


122 


27.4% 


DZ os 


385 


176 


45.7% 


400 


102 


25.5% 


Total 


1505 


595 


39,5% 


1573 


352 


22.4% 
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Table 2 Tetrachoric twin correlations (r) with 95% 
confidence interval (95%CI) for 1 8-25 year old twins from 
two birth cohorts 





1993-1995 


2009-2010 


Combined 
(1993-1995 & 
2009-2010) 




r 


95%CI 


r 


95%CI 


r 95%CI 


MZM 


.82 


.70-.90 


.75 


.43 - .91 


.80 .70 - .88 


DZM 


.51 


.30-67 


.37 


.04 - .72 


.47 .29 - .63 


MZF 


.77 


.66-85 


.80 


.65 - .89 


.78 .70 - .85 


DZF 


.69 


.55 -.81 


.66 


.44 - .81 


.67 .55 - .77 


DOS 


.40 


.26 - .54 


.34 


.11 - .54 


.38 .25 - .49 



over the 2 cohorts did not deteriorate the fit of the model 
(Table 4, model 3 for men, model 4 for women, model 5 
for men and women simultaneously), indicating that the 
estimates were not significantly different between the two 
cohorts. Based on the combined data from both cohorts 
55% (95%CI: 36% - 75%) of the total variance in liability to 
smoking is explained by genes, 23% (95%CI: 7% - 39%) by 
shared environmental factors and the remaining variance 
of 21% (95%CI: 16% - 28%) is explained by unique envir- 
onmental factors. 

There was a trend (see Table 3) for heritability to be 
higher in men (62%-67%) than in women (15%-34%) 
while shared environmental influences were higher in 
women (.46-. 62%) compared to men (.06%-19%). The 
test for sex differences in the contribution of genetic 
and environmental factors showed a p-value of 0.0195 
(Table 4, model 6). With a critical p-value of 0.01 we 
concluded that the sex differences were not significant. 

Table 3 Genetic architecture of smoking behavior 



1993-1995 2009-2010 Combined 

(1993-1995 & 
2009-2010) 





Estimate 


95% 
CI 


Estimate 


95% 
CI 


Estimate 


95% 
CI 


Man: 


Genetic factors 


.62 


.23- 
.86 


.67 


.00 - 

.90 


.66 


.29 - 
.86 


Shared 
environment 


.19 


.01- 
.53 


.06 


.00 - 
.66 


.14 


.00 - 
.46 


Unique 
environment 


.18 


.10- 
.30 


.26 


.09 - 
.57 


.20 


.12 - 

.30 


Woman: 


Genetic factors 


.15 


.00- 
.49 


.34 


.00 - 
.84 


.21 


.00 - 

.49 


Shared 
environment 


.62 


.32- 
.80 


.46 


.00 - 

.79 


.57 


.32 - 
.77 


Unique 
environment 


.23 


.15- 
.33 


.20 


.10 - 
.34 


.22 


.15 - 

.30 



Estimates for the contribution of additive genetic, common environmental 
and unique environmental factors to the variance in smoking in 18-25 year 
old twins in 1993-1995 and 2009-2010 and in the combined sample (full 
model) 



Dropping the genetic factors or the shared environmen- 
tal factors from the model caused a significant deteriora- 
tion (Table 4, respectively model 7 and 8) in goodness- 
of-fit. 

Discussion 

This study explored the prevalence and heritability of 
smoking initiation in two birth cohorts of young adults 
(18-25 years) and tested for Genotype x Environment 
interaction by estimating heritability conditional on 
birth cohort. The first cohort was born in 1968-1977 
and their smoking behavior was assessed in 1993-1995. 
The prevalence of smoking was 50.9% in men and 
39.5% in women. The second cohort was born in 1984- 
1992 and participated in the longitudinal survey study of 
the Netherlands Twin Register in 2009-2010. In this 
cohort a substantially lower prevalence of smoking was 
seen in both sexes with a prevalence of 22.8% in men 
and 22.4% in women. Thus, significantly more young 
adults started smoking in 1993-1995 compared to sub- 
jects of the same age in 2009-2010. Those results prob- 
ably reflect the change in smoking policy in the 
Netherlands, indicating that the discouragement of 
smoking is successful. The observed decrease is in line 
with other studies in the Netherlands (a decrease of 
smoking from 38% to 30% for man and from 30% to 
26% for women, from 1991 to 2009, http://www.stivoro. 
nl) suggesting that the NTR sample is representative for 
the Dutch population. A decrease in the prevalence of 
smoking is also seen in other European countries (a 
decrease in prevalence in Swiss and Finnish men [2,3]) 
and in the United States (a decline of smoking initiation 
among high school students [1]). 

We evaluated the genetic architecture of smoking in 
both cohorts and did not find changes in the heritability 
as a function of environmental exposure. There are dif- 
ferent theories to explain GxE interaction mechanisms. 
Policies and rules with respect to smoking behavior 
have changed between 1993-1995 and 2009-2010. This 
changed the general attitude from smoking as acceptable 
behavior to antismoking (= social push), and could 
increase heritability since only genetic vulnerable sub- 
jects still smoke. On the other hand, when smoking is 
stigmatized and banned in public places, social forces 
affect smoking among most persons, both with and 
without genetic tendencies to smoke (= social control). 
These social controls will cause an actual drop in smok- 
ing among genetically motivated smokers and will 
decrease the genetic influence on smoking [13]. The fact 
that we did not observe any changes in the heritability 
of smoking between both cohorts suggests either that 
those two mechanisms (social control and social push) 
counterbalance each other or that there are no GxE 
interactions for smoking in this age group. 
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Table 4 Comparison of ACE models in both birth cohorts 


Compare cohorts 


-2 LL 


Df 


Versus 


X 2 


a df 


P 


I Full ACE 


5676.915 


4990 










2 Constrain Rcdos at 1 


5677.127 


4992 


I 


0.212 


2 


0.8994 


3 EQ ACE for men over cohorts 


5677.825 


4995 


2 


0.698 


3 


0.8737 


4 EQ ACE for women over cohorts 


5677.692 


4995 


2 


0.565 


3 


0.9044 


5 EQ ACE for both sexes over cohorts 


5678.229 


4998 


2 


1.102 


6 


0.9814 


6 No sex differences in contributions of A,C,E 


5688.1 1 7 


5001 


5 


9.888 


3 


0.0195 


7 Drop A 


5718.180 


5002 


6 


39.951 


1 


< 0.0001 


8 Drop C 


5695.529 


5002 


6 


7,412 


1 


0.0064 



Full ACE: the influence of A, C and E factors is estimated separately for men and women and for both cohorts. The shared environmental correlation (Rc) in DZ 
twins of opposite sex (DOS) is estimated as a free parameter. Model 2: test for qualitative sex differences (Rcdos is fixed at 1). Model 3: ACE in cohort 1993-1995 
is equalized to ACE in cohort 2009-2010 for men. Model 4: ACE in cohort 1993-1995 is equalized to ACE in cohort 2009-2010 for women. Model 5: ACE in cohort 
1993-1995 is equalized to ACE in cohort 2009-2010 for men and women simultaneously. Model 6: ACE for men is equalized to ACE for women. Model 7: Additive 
genetics effects (A) are dropped. Model 8: Shared environmental influences (C) are dropped 



In a Swedish study, regular tobacco use was explored in 
3 different birth cohorts: 1910-1924, 1925-1939, 1940- 
1958. For men, the prevalence (65%, 66%, 72%) and the 
heritability estimates (53%, 58%, 51%) all were stable 
across the 3 birth cohorts. In contrast, the prevalence in 
women increased over time, from 18% in 1920-1924 to 
52% in 1940-1958. The heritability estimates also increased 
remarkably from 0% to 64%. Kendler et al. suggested that 
heritability increased in women because the social restric- 
tions on tobacco use in women became more relaxed over 
time [17]. The birth cohorts in the Swedish study were 
much earlier than the cohorts participating in the NTR 
study. The different outcomes in men and women in the 
Swedish data are therefore not surprising. 

Heath et al (1993) explored data on smoking initiation 
in three adult twin samples from Australia, Virginia and 
AARP. Data were broken down into 4 or 5 birth cohorts 
ranging from 1915-1939 to 1960-1967. Despite a marked 
change in the proportion of male (not female) respon- 
dents who reported ever having smoked, no evidence 
was found for cohort differences in genetic and environ- 
mental effects (no Genotype x Cohort interaction) for 
either men or women [24]. Those results are in line 
with the present study. Heath et al. suggest that the 
same variables, like personality traits of biological 
effects, have continued to determine which individuals 
will become smokers, but that in more recent years only 
individuals with more extreme values on those traits are 
becoming smokers. 

In the studies of Kendler et al and Heath et al, birth 
cohort and age effects are intertwined. Kendler et al 
showed that 20% of the total variance in lifetime regular 
tobacco use was explained by age. In contrast, in our 
study we analysed data of two groups of twins from dif- 
ferent birth cohorts both measured when they were 
between 18 and 25 years old. 

A study from the USA explored the heritability for 
smoking onset in different states rather than for 



different birth cohorts (mean age 16.4 years, age range 
12-21, n = 2060 twin- and sibling pairs from 31 states) 
and reported that the average rate of smoking onset var- 
ied significantly across states (overall, 45% of the 
respondents reported smoking onset) but the measure 
of genetic influence did not. In contrast, the genetic 
influence for daily smoking varied across states [25]. 

Although the samples described in the other studies 
have different backgrounds compared to ours, the 
results suggest that the heritability of smoking initiation 
is not influenced by changes in prevalence ([24,25] and 
present study) while for regular/ daily smoking the herit- 
ability seems to vary with fluctuating prevalence [17,25]. 

Boardman suggested that genes for smoking initiation 
may be different from those associated with daily smok- 
ing. Initiation of smoking might be more heavily influ- 
enced by personality characteristics such as novelty 
seeking or impulsivity, whereas regular tobacco use is 
related to other mechanisms such as nicotine metabo- 
lism. In line with this theory, recent meta-analyses of 
genome-wide association studies for smoking pheno- 
types detected significant associations between number 
of cigarettes smoked per day and polymorphisms in 
nicotine metabolism genes while those associations were 
not detected for smoking initiation [26,27]. 

We found that the genetic architecture for smoking 
initiation did not differ as a function of environmental 
conditions. In the data of the NTR, we also observed a 
change in the prevalence of alcohol use over the past 15 
years (an increase). The heritability for alcohol use 
(alcohol initiation, frequency and quantity) was stable 
over the cohorts (Geels LM, Bartels M, van Beijsterveldt 
CEM, van der Aa N, Vink JM, Boomsma DI: Trends in 
adolescent alcohol use: Effects of age, sex and cohort, 
submitted). Those results are in line with our results on 
smoking initiation and suggest that changes in environ- 
mental factors causing a change in prevalence not auto- 
matically lead to change in the genetic architecture. 
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Our full genetic model suggested higher heritability 
estimates for men than for women while shared envir- 
onmental influences were higher in women than men. A 
formal test showed that a model without sex differences 
did not worsen the fit of the model significantly, based 
on a p-value threshold of 0.01 but this test was border- 
line significant (p = 0.0195). Other studies exploring sex 
differences in genetic architecture showed contrasting 
results. A Finnish study on smoking initiation (early, 
late, never) reports a higher influence of genetic factors 
in men compared to women [7] while a meta-analysis 
for smoking initiation (based on 6 studies) showed a 
higher heritability in women compared to men [11], 
suggesting that if sex differences in heritability exist, 
they probably are not large. We did not find sex-differ- 
ences in heritability of smoking initiation in an adult 
NTR sample [4] with a mean age of 30.5 years. The 
cohorts in the present study consist of younger adults 
(18-25) and at this age the contribution of environmen- 
tal influences could differ for men and women. How- 
ever, in a previous study of twins registered with the 
NTR, the predictors of the uptake of smoking were 
similar in boys and girls [28]. 

Conclusions 

In summary, these results describe a remarkable decrease 
in smoking between 1993-1995 and 2009-2010, but no 
changes in the relative contribution of genetic and envir- 
onmental factors. Variation in smoking in 18-25 year 
olds was explained by genetic (55%), shared environmen- 
tal (23%) and unique environmental (21%) factors. Those 
results are in line with other twin studies, both in Dutch 
samples and worldwide [4,5,9,11,29-31]. Future studies 
should incorporate additional phenotypes to test whether 
heritability will change for other smoking related pheno- 
types, such as number of cigarettes smoked per day or 
nicotine dependence. Furthermore, this research should 
be extended to other age groups. We focused on young 
adult twins aged 18-25 years because adolescence and 
young adulthood is the critical period to start smoking 
and relatively few individuals will start smoking after the 
age of 25. Finally, future studies of smoking behavior 
including genome-wide association studies should not 
only include age and sex, but should also consider birth- 
year and/or year of assessment as important covariates to 
explain differences in prevalence. 
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